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I.  INTRODUCTION 


The  focus  of  this  research  has  been  primarily  on  the  development 
and  testing  of  interactive  methodology  for  attacking  very  complex  and 
difficult  scheduling  and  routing  problems.  The  research  has  involved 
three  major  components;  (1)  determining  the  level  of  human  interaction, 

(2)  developing  an  interactive  interface,  and  (3)  developing  mathematical 
models  to  aid  in  the  decision  making  process.  The  interactive  interface 
has  been  built  around  a  colorgraphics  computer  terminal.  Two  such  units 
have  been  purchased  by  the  School  of  Industrial  and  Systems  Engineering 
at  Georgia  Tech  as  a  result  of  research  sponsored  under  this  contract. 
Developing  the  interactive  interface  and  determining  the  level  of  human 
interaction  have  been  the  most  time  consuming  activities  since  there  has 
been  very  little  previous  research  on  which  to  build  in  this  area. 

While  we  have  found  that  to  be  effective,  interactive  methoao^>gy 
must  be  tailored  to  the  particular  application,  there  are  some  unifying 
underlying  concepts.  The  first  is  that  some  spatial  representation  of  the 
problem  is  essential  if  the  human  is  to  play  a  major  role  in  the  optimization 
phase  of  the  decision  making  process.  If  the  role  of  the  human  only  involves 
selecting  from  a  few  alternatives  generated  by  an  optimization  model  or 
heuristic,  numerical  information  may  suffice.  However,  if  the  human  is 
to  help  guide  the  optimization  process,  human  insight  regarding  the  optimi¬ 
zation  process  is  much  better  when  aided  by  some  spatial  representation. 

For  scheduling  and  routing  problems,  we  have  found  bar  charts  of  resources 
versus  time  and  network  representations  to  be  very  insightful.  While  these 
representations  include  numerical  information  (e.g.,  time,  capacity,  distance, 
ets.)  the  primary  insight  seems  to  come  from  the  picture  rather  than  the 
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numbers . 

In  developing  our  methodology,  we  have  looked  at  three  specific  problem 
areas:  (1)  long  range  fleet  scheduling,  (2)  vehicle  delivery,  and  (3)  order 
picking.  These  topics  are  discussed  in  detail  in  Sections  II,  III,  and  IV. 
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II.  FLEET  SCHEDULING 

The  long  range  fleet  scheduling  problem  is  a  very  important  but  very 
difficult  problem.  It  is  difficult  from  a  mathematical  point  of  view  be¬ 
cause  it  has  a  very  complex  constraint  structure  and  a  huge  number  of 
possible  schedules.  It  is  also  difficult  from  a  philosophical  point  of 
view  since  the  objective  is  not  precisely  defined  and  tends  to  varv  as  con¬ 
ditions  change.  The  fleet  scheduling  problem  can  best  be  described  in 
terms  of  a  specific  class  of  ships.  We  will  talk  in  terms  of  the  cruisers 
and  destroyers  which  are  a  part  of  the  Atlantic  Fleet. 

Very  briefly  the  SURFLANT  long-range  CRU/DES  scheduling  problem 
is  as  follows.  Roughly  70  destroyers  and  10  cruisers  are  to  be  scheduled 
over  a  60  month  planning  horizon.  There  are  a  number  of  "hard"  commitments 
(e.g.,  SIXTHFLT,  NATO,  etc.)  which  must  be  satisfied.  There  are  other 
commitments  (e.g.,  service  to  other  fleet  training)  which  are  "soft"  in 
the  sense  that  they  are  satisfied  when  possible  but  not  at  the  expense 
of  the  hard  commitments.  Each  unit  has  a  periodically  scheduled  overhaul. 
There  are  numerous  "constraints"  on  deployments.  For  example,  a  unit 
should  not  be  redeployed  for  at  least  six  months  following  a  deployment 
with  eight  or  nine  months  desirable.  Units  require  a  workup  period- of 
six  or  seven  months  after  overhaul  before  deployment  with  eight  or  nine 
months  desirable.  Deployment  groups  should,  if  possible,  have  geographical 
integrity.  While  the  above  properties  are  incomplete,  they  do  indicate 
the  flavor  of  the  scheduling  problem. 

This  problem  is  currently  being  solved  using  pencil  and  paper  methods. 
There  are  some  obvious  drawbacks  to  using  manual  methods  to  attack  scheduling 
problems  of  this  size  and  complexity.  Some  of  these  shortcomings  are  vividlv 
indicated  by  LCDR  Bobst's  comments  in  Appendix  C  of  [3].  Manual  methods 
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are  very  time  consuming,  hence,  it  is  difficult  to  react  to  changes  in  the 
system.  A  great  deal  of  insight  into  the  scheduling  problem  is  required. 
This  insight  is  very  difficult  to  document  and  in  many  cases  can  be  ob¬ 
tained  only  by  extensive  hands  on  experience.  It  is  impossible  to  consider 
more  than  a  very  small  fraction  of  the  possible  schedules.  And  perhaps 
most  importantly,  in  cases  where  the  system  constraints  are  "tight"  a 
scheduler  may  not  be  able  to  generate  a  feasible  schedule-  even  after  many 
hours  of  effort. 

Given  these  difficulties  with  manual  scheduling  methods,  it  seems 
natural  to  consider  computerization  of  at  least  a  portion  of  the  scheduling 
function.  The  first  level  of  computerization  is  obviously  some  form  of 
"information  system."  This  is  in  simple  terms  a  bookkeeping  system  for 
keeping  track  of  the  various  parameters  relevant  to  the  scheduling  problem. 
One  such  system,  which  is  apparently  operational,  is  described  in  DPSCPAC 
documents  [1J  and  [2],  The  report  writing  capability  of  this  system  allows 
information  about  the  schedule  to  be  printed  out  in  various  forms  including 
a  sort  of  modified  Gantt  chart.  The  system  is  not,  however,  a  schedule 
generator.  The  schedule  must  be  inputted  by  the  user.  The  system  has  the 
flexibility  to  allow  the  user  to  create  schedules  from  a  set  of  standard 
modules  but  the  modules  themselves  must  be  provided  by  the  user.  The  system 
being  proposed  in  OKI  report  [3]  is  also  an  information  system  rather  than 
a  schedule  generator. 

These  systems  have  some  advantages  over  manual  methods.  In  particular, 
they  allow  the  user  to  very  rapidly  access  information  about  given  schedules 
and  they  allow  this  information  to  be  organized  in  a  variety  of  forms. 

They  force  a  degree  of  "standardization"  which  may  be  desirable.  And,  in 


theory  at  least,  they  make  it  possible  for  the  scheduler  to  more  easily 
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check  schedule  feasibility  and  hence  allow  the  consideration  of  more 
potentially  good  schedules. 

They  have  as  potential  disadvantages  the  requirements  that  the  user 
must  be  knowledgeable  about  the  computer  system  in  addition  to  being  know¬ 
ledgeable  about  the  scheduling  problem  itself.  Also,  the  user  must  obviously 
have  access  to  a  computer.  In  order  to  have  any  real  effectiveness,  this 
access  needs  to  include  an  interactive  capability.  To  some  extent,  these 
difficulties  will  be  encountered  in  any  form  of  computer  aided  scheduling. 

The  greatest  shortcoming  of  information  systems  of  the  type  indicated 
above  is  that  they  only  have  the  capability  to  speed  up  the  scheduling 
process.  They  do  not  bring  to  bear  any  of  the  available  mathematical 
power  to  overcome  problems  the  user  must  face  in  resolving  the  fundamental 
combinatorial  structure  inherent  in  the  scheduling  problem. 

Mathematical  Model 

In  order  to  determine  the  respective  roles  of  the  computer  model  and 
the  human  in  an  interactive  process,  it  must  first  be  ascertained  how  much 
of  the  problem  structure  can  be  captured  in  a  computationally  tractable 
mathematical  model.  After  examining  a  number  of  alternative  formulations, 
it  was  determined  that  a  generalized  set  partitioning  model  was  the  most 
appropriate  for  modeling  this  problem. 

In  order  to  illustrate  the  generalized  set  partitioning  structure 
consider  the  following  simple  example.  (The  structure  generalize  in  a 
straight  forward  fashion  as  will  be  indicated  later.)  Suppose  that  there 
are  five  destroyers  to  be  scheduled.  Destroyers  1  and  2  are  based  at 
Charleston  while  destroyers  3,  A,  and  5  are  based  at  Norfolk.  Suppose  that 
there  is  a  sjngle  commitment  to  be  satisfied.  Suppose  that  the  commitment 
is  to  provide  two  destroyers  at  all  times.  Let  the  planning  horizon  be 
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twelve  months. 

We  will  enumerate  the  possibilities  for  each  destroyer  in  a  sort  of 
expanded  Gantt  Chart  as  illustrated  in  Figure  1.  Suppose  that  Destroyer  1 
can  complete  all  of  its  desired  workup  and  then  deploy  during  April  through 
September.  (The  possible  deployments  are  represented  in  rows  1  through 
12  of  Figure  1.)  '[’his  deployment  is  represented  by  column  x^.  Another 
possibility  is  to  decrease  the  desired  workup  by  one  month  and  deploy 
destroyer  1  during  March  through  August.  This  is  represented  by  colunm 
X12‘  ^  third  possiblity  is  to  decrease  the  desired  workup  by  two  months 

and  deploy  destroyer  1  during  February  through  July.  This  is  represented 
by  column  x^-  Constraint  13  insures  that  we  don't  choose  more  than  one 
of  the  possiblities  for  destroyer  1. 

Suppose  that  destroyers  3  and  4  are  currently  deployed.  They  can  con¬ 
tinue  their  deployment  through  February,  have  six  months  off  and  be  rede¬ 
ployed  during  September  through  December.  These  deployments  are  represented 
by  column  x^j  for  destroyer  3  and  x^  for  destroyer  4.  Another  alternative 
might  be  to  cut  a  month  off  their  current  deployment  (i.e.  continue  their 
deployments  through  January),  give  them  six  months  off  and  redeploy  them 
during  August  through  December.  These  deployments  are  represented  by 
columns  x ^  and  x^?.  Similarly,  deployments  for  destroyers  2  and  5  are 
represented  by  columns  x^,  x22’  X51’  X52"  anc^  X53’ 

By  utilizing  this  structure,  we  can  include  most  of  the  essential 
elements  of  the  scheduling  problem.  Having  a  longer  planning  horizon  and 
more  than  one  requirements  area  simple  increases  the  size  of  the  problem. 
There  are  obviously  some  complex  trade-offs  which  may  have  to  be  made  in 
the  schedule  generation  process.  For  example,  it  may  not  be  possible  to 
generate  a  schedule  which  has  both  the  property  that  deployment  groups 
have  area  integrity  and  each  unit  has  the  desired  workup  period  prior  to 
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Figure  1:  Illustration  of  the  Integer  Program 
Underlying  the  Scheduling  Generator 
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deployment.  To  handle  such  trade-offs  in  the  model,  a  system  is  required 
for  weighting  the  alternatives.  These  weights  can  then  be  used  to  "drive" 
the  schedule  generator. 

Since  there  is  a  column  in  t fie  model  corresponding  to  every  possible 
individual  ship  schedule,  there  are  potentially  a  very  large  number  of 
columns.  The  tractability  of  the  model  is  primarily  a  function  of  the 
number  of  columns,  lienee,  it  was  necessary  to  determine  the  approximate 
number  of  columns  one  could  expect  in  an  application.  Unfortunately,  this 
turned  out  to  be  a  very  difficult  task.  The  only  way  to  accurately  make 
this  determination  is  to  actually  generate  the  possible  schedules  using 
realistic  data. 

Some  historical  data  was  obtained  with  the  help  of  the  SURFLANT  CRU/DES 
schedulers.  Using  this  data  a  schedule  generator  was  developed  and  tested 
to  determine  the  number  of  potential  possible  schedules.  The  generator 
included  realistic  features  such  as  the  ability  to  not  allow  consecutive 
deployment  to  certain  area  such  as  the  Middle  East,  minimize  and  maximum 
times  between  deployments,  and  different  overhaul  cycles.  The  sample 
historical  data  included  57  possible  deployments  over  a  four  year  planning 
horizon.  The  planning  horizon  was  divided  into  96  two-week  periods. 

In  general  we  found  that,  as  expected,  the  number  of  possible  schedules 
for  each  ship  varied  depending  on  the  minimum  and  maximum  time  allowed 
between  deployments  and  the  timing  of  the  overhaul  during  the  period. 

For  the  data  we  examined  the  ships  each  had  one  overhaul  during  a  four 
year  period.  It  it  occurred  at  the  beginning  or  end  of  the  four  year 
planning  horizon  then,  there  were  more  possible  schedules  than  if  it  occurred 
in  the  middle.  Assuming  that  the  time  between  deployment  was  between 
five  and  ten  months,  the  average  number  of  schedules  for  the  cases  con¬ 
sidered  was  139  with  a  high  of  165  and  a  low  of  92.  This  means  that  in 
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scheduling  80  ships,  if  all  possible  schedules  for  each  individual  ship 
were  generated,  the  model  would  have  about  10,000  columns. 

Since  this  is  quite  large  for  an  integer  program,  even  one  with 
the  special  set  partitioning  structure,  we  felt  it  necessary  to  look  at 
heuristic  as  welL  as  optimum  methods  for  solving  the  model.  These  methods 
are  based  on  a  pricing  theorem  and  are  essentially  the  same  for  both 
scheduling  and  routing  problems.  This  will  be  discussed  in  detail  in 
Section  III. 

In  general,  we  found  the  methods  that  we  developed  for  solving  the 
set  partitioning  problem  to  be  very  satisfactory.  They  provided  high 
quality  answers  in  a  very  reasonable  amount  of  time.  (liven  the  imprecise 
nature  of  the  values  or  weights  assigned  to  each  column,  we  feel  that 
the  effort  required  to  allow  the  methods  to  verify  optimality  is  not 
justified. 

Weighting  Altern ative  Schedules 

In  utilizing  the  set  partitioning  model,  one  of  the  essential  elements 
is  a  scheme  for  determining  the  value  or  "weight"  which  should  be  given 
to  each  column.  In  fact,  any  scheme  which  could  be  done  without  human 
interaction  would  require  such  a  weighting  scheme.  The  SURFLANT  schedulers 
who  have  had  one  original  contact  at  Norfolk  felt  that  such  a  scheme  could 
be  developed.  However,  as  we  attempted  to  develop  such  a  scheme  and 
discussed  the  concept  with  CMDR  Lee  Pittman,  who  had  recently  taken 
over  the  CRU/DES  scheduling  activity  at  SURFLANT,  it  became  clear  that  a 
general  scheme  was  not  possible.  This  was  due  to  both  the  continuously 
changing  nature  of  the  problem  and  the  level  of  judgement  required  to 
evaluate  each  alternative  schedule.  We  determined  that  the  onlv  realistic 
way  to  evaluate  an  individual  ship  schedule  was  to  have  an  experienced 
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scheduler  actually  look  at  it.  To  generate  all  of  the  10,000  possible 
individual  schedules  and  have  a  person  evaluate  each  one  was  obviously 
not  practical.  This  lead  us  to  look  at  interactive  schemes  where  a 
person  would  help  in  the  generation  of  individual  schedules  (i.e.,  columns 
in  the  model).  By  performing  this  column  generation  interactively,  some¬ 
thing  is  lost  in  terms  of  optimality  since  only  a  small  subset  of  the 
possible  columns  would  ever  be  generated.  However,  the  person  will 
hopefully  be  able  to  generate  very  good  candidate  columns  which  will  then 
result  in  a  good  overall  schedule. 

Since  the  underlying  problem  structure  and  the  set  partitioning 
model  are  essentially  the  same  for  both  the  fleet  scheduling  problem 
and  a  broad  class  of  delivery  problems,  it  was  decided  to  address  both 
problem  classes  in  developing  the  interactive  methodology.  This  had  two 
distinct  advantages.  There  were  results  available  in  the  literature  on 
the  delivery  problem  against  which  an  interactive  approach  could  be  compared 
while  none  were  available  for  the  fleet  scheduling  problem.  Also,  the 
delivery  problem  is  generally  simpler  than  the  fleet  scheduling  problem. 

We  felt  that  some  experience  was  needed  in  developing  interactive  methodology 
for  the  simpler  problem  before  attacking  the  more  complex  fleet  scheduling 
problem. 

The  methodology  developed  for  the  delivery  problem  is  discussed  in 
Section  III.  We  feel  that  this  represents  a  real  breakthrough  in  attacking 
this  class  of  problems.  A  further  refinement  of  this  methodology  is  being 
continued  under  a  current  ONR  contract  and  will  be  reported  on  in  the  near 
future.  We  felt  that  this  same  basic  methodology  with  a  different  inter¬ 
active  interface  could  be  extended  to  solve  the  fleet  scheduling  problem. 
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FI  eet  Scheduling  Interactive  Interface 

The  first  interface  that  we  considered  was  a  bar  chart  color  coded 
to  indicate  the  various  deployment  areas.  Each  ship  is  represented  as  a 
bar  with  time  periods  (two  weeks)  along  the  horizontal  axis.  One  example 
of  the  colorgraphics  display  is  shown  in  Figure  2.  This  system  allows  the 
user  to  add  or  delete  deployments  to  a  given  ship’s  schedule,  to  move  de¬ 
ployments  from  one  ship  to  another,  to  specify  overhaul  cycles,  and  to 
generate  various  statistics  such  as  minimum  or  average  time  between  de¬ 
ployment  for  a  given  ship.  This  is  very  similar  to  the  manual  bar  charts 
currently  used  to  generate  schedules.  While  certain  features  such  as  the 
abilitv  to  generate  statistics  were  very  helpful,  overall  the  display  did 
not  convey  enough  information  to  allow  the  human  to  bring  his  insight  to 
bear  on  the  problem.  In  particular  it  was  designed  to  handle  individual 
ships  while  the  scheduler  is  most  frequently  concerned  with  developing 
battle  groups. 

Hence  a  second  and  much  more  elaborate  interface  was  designed.  We 
proposed  that  the  entire  scheduling  function  be  translated  into  a  mini¬ 
computer-based  interactive  scheduling  program  package  and  that  this  package 
be  designed  to  use  the  extensive  colorgraphics  capabilities  of  the  Chromatics 
CG  1999  minicomputer.  The  package  would  consist  of  four  major  programs 
and  a  master  module  with  menu  selection  from  the  master  to  the  four  major 
proqrams.  The  modules  are  described  in  Figure  3. 

Ship  File  Manager.  This  module  is  used  to  maintain  a  file  on  diskette 
of  all  hulls  in  the  Fleetpak  data  base.  The  program  functions  supported 
are  List,  Displav,  Add,  Delete,  Change,  and  Return.  Table  1  shows  the 
list  of  data  elements  that  will  be  maintained  for  eacti  ship.  Note  that 
the  initial  version  of  Fleetpak  contains  only  the  first  seven  data  elements. 
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-SFMGR  ship  file  manager 
-DFMGR  deployment  file  manager 
-BGMGR  battle  group  manager 
-PSMGR  permanent  schedule  manager 


-List 

-Display 

-Add 

-Delete 

-Change 

-Return 


-List 

-Display 

-Add 

-Delete 

-Change - Change  Name 

-Return  Change  Dates 

Add  Slot 
Delete  Slot 
Change  Slot 
Return 


-Save 

-Load 

-Alter - Search 

-Return  Relax 

Restrict 

Assign 

Release 

Return 


-List - Ship  Schedule 

-Display - Ship  Data 


-Adopt  Deployment 
-Drop  Deployment 
-Adopt  Battle  Group 
-Drop  Battle  Group 
-Return 


Figure  J.  Fleetpak  Structure 
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Table  1.  Ship  File  Contents 


Data  Element 

Bytes  Required 

Ship  Tupe 

4 

these  elements 

Hull  Number 

4 

in  initial  version 
of  Fleetpak 

Ship  Name 

25 

Ship  Class 

25 

Characteristic  Word 

U1C 

Home  Port  CO 

Home  Port 

CO 

ADSCN 

8 

Displacement 

Length 

Beam 

Draft 

Officers 

Enlisted 

Sq  Ft  Capacity 

Cu  Ft  Capacity 

Max  Speed 

Econ  Speed 

Max  Speed  Range 

Econ  Speed  Range 

Propulsion 

Number  of  Shafts 

Fuel  Type  1  Fuel  Type  1  Capacity 

Fuel  Type  2  Fuel  Type  2  Capacity 

Fuel  Type  3  Fuel  Type  3  Capacity 

Fuel  Type  4  Fuel  T^ge_4_Cagacitj;_ 

Number  on  Board  _  _ for_all_sys terns 
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Later  versions  may  support  the  entire  list. 

The  Fleetpak  user  would  necessarily  use  the  Ship  Kile  Manager  as  the 
first  step  in  building  a  schedule.  All  the  information  about  available 
ships  would  be  entered.  Later,  changes,  additions,  and  deletions  would 
be  processed.  At  any  time,  information  about  an  individual  ship  can  be 
displayed  or  listed.  A  list  of  all  ships  would  also  be  provided. 

The  logical  structure  of  the  ship  file  is  a  stack.  The  physical 
organization  is  sequential.  Access  methods  supported  are  sequential  and 
hashed  lookup  with  sequential  lookup  of  overflows  in  an  overflow  table 
at  the  end  of  the  file. 

Deployment  File  Manager.  This  module  is  used  to  maintain  a  file  on 
diskette  of  all  deployments  in  the  data  base.  The  program  supports  List, 
Display,  Add,  Delete,  Change,  and  Return  functions.  Changes  possible 
include  Change  Name,  Change  Dates,  Add  Slot,  Delete  Slot,  Change  Slot, 
and  return.  The  Change  Slot  function  includes  Change  Type,  Change,  Class, 
Add  Characteristic,  and  Delete  Characteristic. 

The  deployment  file  manager  program  is  used  to  define  all  the  available 
ship  activities  to  which  ships  could  be  assigned  singly  or  in  battle  groups. 

Table  2  shows  the  contents  of  the  Deployment  File.  Table  3  show  the 
related  file  called  tie  Event  File.  No  direct  user  access  to  the  Event  File 
is  possible.'.  It  is  an  incex  file  used  by  the  various  modules  of  Fleetpak. 
The  Deployment  File  provides  access  to  schedule  data  by  deployment  name. 

The  Event  File  provides  access  bv  ship  type:  hull  number.  Two  files  are 
needed  to  give  quick  access.  lor  example,  it  would  be  very  slow  to  try 
to  find  all  the  deployments  scheduled  for  an  individual  ship  by  reading 
through  all  deployments  and  looking  within  them  for  an  entry  for  the  in¬ 
dividual  ship. 


r 
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Table  2.  Deployment  File  Content;; 


Data  Element 
Deployment  Name 
Description 
Start  Year 
Start  Month 
Start  Day 
End  Year 
End  Month 
End  Day 

Long  Range/Short  Range  Code 
LR/SR 

Number  of  Slots 
Ship  Type 
Ship  Class 

Required  Character istics  Word 


Bytes  Required 
8 
40 
2 
2 
2 
2 
2 
2 
6 
1 
2 
4 

1  entry  for 

^  each  slot 

g  - 


. .  jv  .«*. .  n  ■wiiimi  •!  ii— i  I,  _  *£,  ,  v  -4r — 

/  -  '  - -  - 
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Table  3.  Event  File  Contents 


Data  Element 


Bytes  Required 


Ship  Type 


Hull  Number 

Number  of  Deployments 

Deployment  Name 

Required  Characteristics  Word 

Start  Year 

Start  Month 

Start  Day 

End  Year 

End  Month 

End  Day 


1  entry  for 


deployment 
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Battle  Group  Manager.  The  Battle  Group  Manager  module  is  the  most 
important  module.  The  purpose  of  BGMGR  is  to  give  interactive  assignment 
of  available  ships  to  open  slots  in  scheduled  deployments.  This  is  the 
third  step  in  building  a  schedule.  An  empty  deployment,  essentially  a 
set  of  slots,  is  called  up  and  reviewed.  For  each  slot,  ships  are  selected 
that  match  the  characteristics  selected  for  that  slot.  In  addition,  the 
ships  must  be  available  during  the  period  of  the  deployment.  To  assist 
the  operator, BGMGR  provides  the  capability  of  displaying  any  individual 
ship  schedule,  any  other  deployment  or  battle  group,  or  the  overall  schedule 
of  all  deployments.  The  ship  file  entry  for  an  individual  ship  can  also 
be  recalled  and  displayed. 

The  purpose  of  having  so  many  different  displays  is  to  give  the  user 
the  opportunity  of  knowing  why  ships  have  already  been  given  other  assign¬ 
ments  and  of  making  intelligent  choices  for  including  ships  in  a  current 
battle  group.  Once  a  slot  is  being  filled,  the  operator  selects  the 
Search  function.  All  ships  that  match  the  characteristics  word  are 
displayed  in  column  2  of  the  Battle  Group  Manager.  From  this  list,  the 
operator  makes  a  selection.  If  he  is  unsure,  then  he  can  display 
the  schedule  for  the  ship  in  question  or  look  at  its  activity  in  other 
deployments . 

Within  this  procedure,  the  user  can  either  relax  or  restrict  the 
characteristics  required  for  each  slot.  Once  a  satisfactory  match  is  made, 
an  assignment  is  made.  Changes  can  be  made  by  releasing  the  ship  already 
assigned  to  a  slot. 

Thus,  the  functions  available  in  the  BGMGR  module  are  Save,  Load, 

Alter,  and  Return.  Under  the  Alter  function  are  subfunctions  of  Search, 
Relax,  Restrict,  Assign,  Release,  and  Return. 
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Note  that  the  initial  version  of  Fleetpak  supports  a  one  to  one 
mapping  of  battle  groups  to  deployments.  Later  versions  may  allow  alter¬ 
nate  schedules  where  deployments  and  battle  groups  meet  on  a  many  to  many 
basis.  In  any  event,  a  complete  schedule  consists  of  a  number  of  adopted 
deployments  and  an  equal  number  of  adopted  battle  groups. 

It  is  within  the  BGMGR  Search  and  Assign  functions  that  the  real 
scheduling  is  actually  performed  under  Fleetpak.  All  other  functions  are 
essentially  just  support  and  report  writing  functions.  It  is  at  this 
point  exactly  that  optimization  can  be  introduced. 

Permanent  Schedule  Manager.  The  PSMGR  module  is  used  to  build  and 
report  schedules.  Functions  supported  are  List,  Display,  Adopt  Deployment, 
Drop  Deployment,  Adopt  Battle  Group,  Drop  Battle  Group,  and  Return.  Under 
the  List  and  Display  functions  there  are  seven  subfunctions.  There  in¬ 
clude  Ship  Schedule,  Ship  Data,  Deployment  Schedule,  Deployment  Data, 

Battle  Group  Data,  Permanent  Schedule,  and  Return. 

File  structure  for  the  permanent  schedule  file  is  fixed  record  length 
sequential  access.  The  file  is  very  small. 

The  main  colorgraphics  display  for  this  system  is  illustrated  in 
Figure  4.  From  this  display  the  user  has  the  ability  to  display  and  alter 
characteristics  required  for  each  ship  in  a  battle  group,  display  future 
requirements  for  ships  with  these  same  characteristics,  generate  candidate 
ships  for  each  battle  group,  display  the  past  and  future  commitments  for 
candidate  ships  for  each  battle  group,  generate  the  battle  groups,  generate 
various  statistics  related  to  the  battle  groups  and  to  individual  ships, 
and  change  individual  deployments  as  the  requirement  or  ship  availabilities 
change . 

The  modules  were  developed  for  testing  on  the  Chromatics  CG  1999  color 
microcomputer.  Our  discussion  with  CMDR  Lee  Pittman,  the  SURFLANT  scheduler 
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for  cruisers  and  destroyers  lead  us  to  conclude  that  this  system  contained 
the  essential  elements  for  efficient  fleet  scheduling.  However,  our 
initial  tests  indicate  that  the  volume  of  ship  related  data  required, 
exceeds  the  capacity  of  the  microcomputer  to  process.  The  design  is 
being  modified  so  that  the  major  data  manipulation  function  will  be  per¬ 
formed  on  a  host  computer.  The  fundamental  interactive  methodology  which 
we  have  developed  seems  adequate  to  address  this  problem.  Further  refine¬ 
ments  are  being  made  under  the  current  ONR  contract. 


r 


DPSCPAC  System  81,  "Long  Range  Lmployment  Schedule  System  Command 
Manual . " 

DPSCPAC  System  83,  "Long  Range  Lmployment  Schedule  User's  Manual." 

Operations  Research  Draft  Final  Report,  "Automated  Ship  Operational 
Scheduling  (AOS)  - Concept  Definition  Phase,"  10  August,  1977. 
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II L.  INTERACTIVE  DELIVERY 

We  will  consider  here  a  set  part itioning  based  appioach  tor 
so  I v i n  ;  a  broad  class  01  routing  problems.  The  approac.h  is  designed 
to  tak  ■  advantage  of  a  high  level  of  human  interaction;  the 
current  implementation  is  interactive  via  a.  colorgraphics  display. 
Howevei ,  many  of  the  concepts  discussed  here  could  he  easily  implemented 
in  an  automatic  system. 

The  rout inp.  |  robler:  which  motivated  much  of  this  work  is  what 
is  called  the  static  pick  up  and  delivery  problem.  This  problem 
will  be  discussed  in  detail  it:  lat.r  sections.  It  is  one  of  the  more 
complex  members  ol  the  i  lass  if  touting  problems  which  are  amenable  to 
the  approach  presented  here.  This  class  also  includes  many  practical 
delivery  problems.  In  order  to  introduce  the  underlying  methodology 
which  ptovides  the  basis  for  the  approach,  consider  a  very 
simple  delivery  example.  Assume  that  a  depot  is  located  at  the 
square  box  labeled  D  in  Figure  1.  From  this  depot  a  single  delivery 
is  to  be  made  to  each  of  the  points  represented  by  numbered  circles. 

The  numbers  on  arcs  connecting  the  circles  represents  the  travel  dis¬ 
tance  between  delivery  points.  Assume  also  that  each  vehicle  (e.g. 

truck)  can  deliver  to  a  maximum  of  two  points  on  a  single  trip.  The 
object  ive  is  to  de. ermine  which  vehicle  should  deliver  to  each  point 
and  t  lie  routing  fo  the  vehicles  which  minimizes  the  total  distance 
t  r.i vt  I  l«.,i . 

Each  column  in  the  matrix  of  Table  1  represents  one  possible  vehicle 
route.  For  example,  column  one  represents  a  vehicle  travelling  from  the 
depot  to  delivery  point  (1)  and  returning.  The  Cj  row  indicates  the 


Table  1.  Matrix  Corresponding  to  Routes  in  the 
Delivery  Example  ol  Figure  1 


routu _ ]_ _ 2  1  4_ 5 6 _ 7  8  9  10  11 _ 12  13  14  15 

8  10  4  4  2  14  10  8  8  10  11  12  6  6  5 

1  00001  111  0  0  0000 

010001  0001  1  1  000 
001000100100110 
000100010010101 
00001  0001  001  01  1 


r 


t 
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distinct-  traveled  for  each  trip.  For  example,  column  six  represents 

.1  vc  ncle  proceed  in;;  i  rom  the  depot  to  delivery  point  (1),  on  to  delivery 

point  (2)  and  from  there,  hack  to  the  depot.  The  value  of  o  is  14, 

6 

the  total  distance  traveled  for  this  trip.  By  enumerating  each  of  the 
possibilities,  as  has  been  done  for  this  matrix,  the  problem  becomes 
one  of  selecting  a  set  of  columns  such  that  every  row  is  represented 
in  exactly  one  column  and  the  sum  of  the  costs  of  the  columns  selected 
is  the  smallest  possible.  This  integer  program  is  called  a  "set  parti¬ 
tioning  model". 

The  set  partitioning  model  was  originally  proposed  for  the  delivery 
problem  by  Balinski  and  Quandt  [2J.  The  model  is  very  powerful  In  the 
sense  that  many  realistic  route  constraints  and  cost  functions  can  be 
handled  easily  in  the  column  enumeration  process.  The  obvious  short¬ 
coming  of  tiie  model  is  that  there  are  typically  a  very  large  number  of 
columns  to  be  enumerated  and  the  resulting  Integer  program  is  very 
large.  The  appieach  presented  here  is  heuristic  in  the  sense  that  we 
generate  only  a  subset  of  the  possible  columns  or  routes  and  in  general 

we  do  not  solve  the  set  partitioning  model  to  optimality. 

The  set  partitioning  model  has  two  very  desirable  features  for 
interactive  optimization.  The  first  is  that  any  route  generated  can  be 
included  as  a  column  in  the  model.  This  allows  the  human  interactor 
t"  utilize  his/her  Intuition  and  spatial  preception  as  well  as  a  wide 
spectrum  of  mathematical  techniques  to  generate  new  routes.  The  second 
feature  is  that,  unlike  more  general  integer  programs,  a  feasible 
solution  to  the  set  partitioning  model  provides  the  basis  for  pricing 
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information  which  can  be  used  to  generate  new  candidate  columns. 

We  will  restrict  the  class  of  routing  problems  considered  here 
to  be  those  for  which  any  subroute  of  a  feasible  route  is  also  a 
feasible  route  with  cost  less  than  or  equal  to  the  cost  of  the  original 
route.  By  imposing  this  restriction,  as  long  as  there  is  at  least  one 
1  in  every  row  of  the  partial  set  partitioning  model  that  we  have 
enumerated,  we  can  easily  generate  feasible  partitions.  The  only 
other  restriction  that  we  put  on  the  class  of  routing  problems  is 
that  we  be  able  to  pose  them  in  a  natural  way  as  set  partitioning 
problems.  However,  it  should  be  noted  that  if  there  is  not  a  nice 
spatial  representation  of  the  routing  problem,  the  human  interactor 
is  much  more  restricted  in  his/her  rontr ibution. 

Set  Partitioning  and  Row  Prices 

The  set  partitioning  problem  can  be  stated  as 


n 

minimi ze 

*  =  l  cjxj 
j-i  J  3 

n 

subject  to 

£  a  x  *  1 
j-i  1J  J 

for  i  =  1,2 

x  =  0  or  1 

J 

for  j  =  1,2 

For  the  delivery  example  in  Table  1,  the  set  partitioning 


model  has  the  second  row  as  the  value  of  the  c^  and  rows  three 

through  seven  as  the  values  of  the  a.^.  The  variables  which  are 

set  to  one  in  a  soiution  to  the  set  partitioning  problem  will  be 

k  k 

called  a  "partition."  We  will  denote  a  partition  as  J  =  (j|x  =  1] 


Balas  and  Padberg  [1]  provide  a  recent  survev  of  results  related  to 
set  partitioning  problems. 


28 


A  fundamental  idea  underlying  much  of  the  work  presented  here 
is  the  concept  of  "row  prices." 

Definition:  P  =  (Pj iPj’ • • • *Pm)  f-s  a  sef  °f  feasible  row  prices 
corresponding  to  the  partition  j1  if 


m  1 


l  P.a  .  =  c  for  j« 
i  =1  1  1 J  J 


It  will  be  useful  to  interpret  the  price  p^  as  an  estimate  of  the 
cost  to  satisfy  constraint  i  using  solution  X1.  For  the  delivery  problem, 

Pj  is  then  an  estimate  of  the  cost  of  satisfying  the  requirement  of 
delivery  point  i  using  the  route  corresponding  to  partition  J1. 

Theorem  1:  Given  a  set  of  feasible  row  prices  (p^.p^, . . . »p^)  corresponding 
to  partition  1  with  value  z\  any  other  partition  J2  has  value 


Proof : 


z2  =  z1  -  l  -  l  (p}a.  -c  ; 

}eJ2  i=l  L  ^  J 


m  ,  ra  i 

l  2  <  I  Piaij~Cj)  =  l  Pi  l  7  3ij 
jeJ  i-1  1  13  J  i*l  1  jer  1J 


Since  is  a  partition,  £  a  =  1  for  each  i  *  1,2,. ...m.  Also, 

j  eJ  2  iJ 

since  (p^ , p^ , . . . , p^)  are  feasible  row  prices  correspond ing  to  X*  we 

m  j 

have  l  Pj  =  Z  .  Hence  the  result  follows. 
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Step  (3):  If  1  ~a  =  1  then  set  a  ,  =  0.  (Note  from  the  assumption 

Jed  ik 

of  section  I  that  any  subroute  of  a  feasible  route  is  also 

a  feasible  route  the  new  column  is  legitimate) 

2  2 

Step  (4):  let  J  =  J  U  f  k. }  (i.e.,  put  column  k  In  the  new  partition) 

Step  (5):  helete  from  N  all  j  for  which  a..  =1  and  a..  =1  for  some 

ik  j  ] 

i  =  1 , 2 ,  .  . . ,  m . 

Step  (b):  If  N  =  0  stop.  Otherwise  go  to  step  (2). 

Note  that  under  the  assumption  that  any  subset  of  a  route  is 

also  a  feasible  route  (discussed  in  section  1)  this  procedure  will 

2 

always  terminate  with  J  as  a  partition  although  not  necessarily  a 

better  partition  than  J^.  In  addition,  computation  to  date  indicates 

that  an  optimum  or  near  optimum  solution  to  the  set  partitioning 

problem  is  determined  very  quickly  by  repeated  application  of  the 

potential  savings  heuristic.  The  heuristic  is  repeated 

until  either  optimality  is  proven  (i.e.,  for  some  X  all 
m  k 

I  piai-  ~  c-  f.  0)  or  until  some  specified  number  of  partitions 
i=l  J  J 
has  been  generated. 

To  illustrate  the  potential  savings  hueristic,  consider  again 

the  delivery  example  depicted  in  Figure  1  and  Table  1.  Suppose 
that  we  select  =  {1,2,..., 3}  as  an  initial  partition.  A  set  of 

feasible  row  prices  is  given  in  Table  2.  (The  question  of  how 

to  generate  "good'1  feasible  row  prices  will  be  addressed  in  the 

next  section.)  Tne  corresponding  potential  savings  J  p.a  ,  -  c. 

£=1  *  ^  ^ 

are  also  given  in  Table  2.  Applying  the  potential  savings  heuristic 
and  breaking  ties  by  selecting  the  column  with  the  lowest  index 


yields  the  new  partition  J  «  {6,13,5}.  The  new  partition  has  a  cost 
2  1 

of  Z  =22  as  compared  to  a  cost  Z  *=  28  for  the  initial  solution. 


I 
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I 


2  f  2 

Using  the  row  prices  P  and  potential  savings  l  p^a^  '  c. 

L=1  J  J 

of  Table  2  and  reapplying  the  potential  savings  heuristic  yields  the 

3  3 

partition  J  =  {8,10,5}  which  has  a  cost  of  Z  =  20.  Again,  from 

Table  2  we  find  that  using  the  row  prices  given  potential  savings 

.  £  0  for  j  =  l,2,...,n.  Hence,  from  Corrollary  1  the 


I  pU  .  -  c  < 

1=1  1  ^  3  j  - 

partition  J  is  optimum. 


Row  Pricing 

k 

For  a  given  partition  X  a  set  of  feasible  row  prices  is  obtained 
by  allocating  the  column  cost  c  for  each  jeJ  among  the  rows  having 
a^j  =  1.  For  the  delivery  example,  this  corresponds  to  allocating 

the  trip  cost  among  the  delivery  points  of  the  trip.  When  a  column 

k  k 

jeJ  dontains  onlv  one  a_  =  1,  the  row  price  is  p^  =  c  .  However, 

when  a  column j  contains  more  than  one  a  =  1,  there  are  an  infinite 
number  of  possible  sets  of  prices.  As  an  example  consider  the  parti¬ 
tion  J-*  =  {8,10,5}  for  the  problem  in  Table  2.  Since  column  5  has 

3 

only  a  =  1  the  value  p  *  2  is  unique.  Column  8  has  both  a,  „  =  1 
j , j  5  l,o 

and  a^  g  =  1,  hence  the  cost  Cg  =  8  could  be  allocated  between  rows 
1  and  4  in  an  infinite  number  of  ways.  Similiarly,  column  10  has 
both  a2  =  1  and  a3  =  1,  hence  c^g  =  10  could  be  allocated 

between  rows  2  and  3  in  an  infinite  number  of  ways.  If  we  allocate 

3  3  3  3 

Cg  as  p^  =  4  and  p^  =  4  and  allocate  c^g  as  =  5  and  P3  =  5,  the 

m  3  3 

resulting  J  p.a.  -  c  do  not  indicate  that  J  is  an  optimum  parti- 
i=l  L3  3 

3 

tion.  Hence  the  set  of  prices  P  given  in  Table  2  are  clearly  better 

3 

since  they  do  indicate  that  J  is  on  optimum  partition. 

Ideally,  we  would  like  a  set  of  prices  which  would  drive  the 
potential  savings  heuristic  toward  an  improving  solution  and  would 


/ 


indicate  optimality  when  no  improving  solution  is  possible  (i.e.. 

We  would  like  the  prices  to  be  analogous  to  dual  variables  in  linear 
programming).  Unfortunately,  it  is  easy  to  construct  cases  for  which 
no  such  prices  exist  (i.e.,  any  problem  for  which  the  integer  and 
continuous  solution  differ).  For  the  delivery  problem  and  the  more 
complex  dial-a-ride  problem  (to  be  discussed  later),  allocating  column 
cost  in  proportion  to  the  cost  of  serving  the  delivery  points  one-at- 
a-time  is  intuitively  appealing  and  seems  to  work  very  well.  As  an 

illustration  consider  again  the  partition  J  in  Table  2.  Column 

has  a,  0  =  1  and  a,  „  =  1.  The  cost  of  serving  delivery  point  1  if 
l,o  4,0 

it  is  the  only  point  in  a  trip  is  Cj  =  1.  The  cost  of  serving  delivery 

point  4  if  it  is  the  only  point  in  a  trip  is  c,(  =  4.  The  prices  for 

3  cicg  2  ^4^8 

rows  1  and  4  were  determined  as  p.  - — — - — •  =  5.3  and  p  =  - — — 

1  C1  C4  4  1  C4 

=  2.7.  The  other  prices  in  Table  2  were  determined  similiarly. 

Column  Generation 

For  large  scheduling  and  routing  problems  it  is  generally  not 
practical  to  generate  all  columns  of  the  corresponding  set  partitioning 
model.  The  remainder  of  this  paper  will  be  concerned  with  using  in¬ 
formation  gleaned  from  one  solution  via  Theorem  1  to  generate  a  new 
and  hopefully  better  solution.  This  is  accomplished  by  either  generating 
new  columns,  adding  them  to  the  current  set  partitioning  model,  and 
then  resolving  the  model  or  by  using  the  information  from  Theorem  1 
directly  to  generate  a  new  solution.  In  the  latter  case  it  is  not 
necessary  to  retain  the  columns  of  the  set  partitioning  model.  How¬ 
ever,  if  the  columns  are  retained,  it  is  possible  to  further  improve 
the  solution  by  periodically  solving  the  set  partitioning  model. 
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We  should  note  that  for  the  class  of  scheduling  and  routing  prob¬ 
lems  being  considered  here,  it  is  very  easy  to  generate  an  initial 
solution,  i'or  our  examples  we  use  the  identity  solution  (e.g.,  in 
the  deliverv  problem  this  is  the  solution  which  has  each  vehicle  making 
a  single  de  ivery)  as  our  initial  solution.  However,  any  feasible 
solution  could  be  used  as  the  initial  solution. 

Clearlv  there  is  a  broad  spectrum  of  possible  approaches  that 
one  migh.  use  to  generate  new  columns  and/or  new  solutions  to  the 
set  partitioning  model.  In  fact,  variations  of  many  of  the  heuristics 
which  have  been  aoplied  to  delivery  problems  can  be  used  very  effectively 
in  conjunction  with  Theorem  1.  In  the  next  section  we  will  discuss 

the  use  of  the  Cl  irke  and  Wright  [1]  savings  procedure  in  conjunction 
with  the  delivery  problem.  In  later  sections  we  will  discuss  more 

complex  clustering  and  chaining  heuristics  as  we  have  applied  them 
to  the  dial-a-ride  problem. 

Clarke  and  Wright  Procedure  with  Pricing 

The  "savings"  heuristic  of  Clarke  and  Wright  [3]  is  the  most 
widely  known  of  tie  heuristics  developed  to  date  for  delivery  problems. 
The  algorithm  proceeds  by  calculating  a  savings  for  each  pair  of  deli¬ 
very  points  i  and  j  defined  as 


s0(i,j)  «  2d0i  +  2d.()  -  =  d0.  +  dj0  -  dtj 


which  is  the  savings  in  mileage  of  supplying  delivery  points  i  and  j 
on  the  same  route  as  opposed  lo  supplying  them  individually  directly 
from  the  depot  (d()^  is  the  distance  from  the  depot  to  deliverv  point  i ) . 


/ 
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Routes  are  then  constructed  either  one-at-a-L ime  or  in  parallel  by 

considering  pairs  of  points  in  order  of  decreasing  savings  and  including 

them  in  the  same  route  if  such  a  route  is  feasible. 

Suppose  that  we  consider  once  more  the  delivery  example  in  Figure  1 

and  the  covering  solution  information  in  Table  2.  Note  that  the  Clarke 

m  1 

and  Wright  savings  values  are  exact Ly  the  values  of  [  p.a  -  c.  in 

i=l  1  J  J 

Table  2  since  d  ^  =  p ,/L  for  i  =  l,2,.,.,m.  Applying  the  C-W  savings 

2 

algorithm  yieLds  the  same  routing  configuration  as  J  =  f6,13,5}.  At 

this  point,  the  C-W  algorithm  would  terminate.  However,  suppose  that 
2 

we  set  d^  =  P^/2  for  i  =  l,2,,..,m.  The  new  savings  are  then  exactly 
2 

t  he  >  p  .  a  . 


i  =  l 


i  ij 


c.  in  Table  2.  Applying  C-W  algorithm  yields  the 
,3 


same  routing  con f i gura t i on  as  J  =  {8,10,5}  which  as  noted  previously  is 
the  optimum  solution  to  this  delivery  example.  Hence,  for  this  example 
at  least  the  C-W  algorithm  without  pricing  did  not  yield  an  optimum  solution 
while  the  same  algorithm  when  combined  with  pricing  did  vield  the  optimum 
solut ion . 

Note  that  wh>  n  routes  are  allowed  to  contain  more  than  two  trips, 
only  a  small  subsol  of  '.he  set  partitioning,  columns  would  be  generated 
using  this  procedure.  \lso,  note  that  the  prices  and  savings  can  be 
calculated  without  ever  generating  the  set  partitioning  matrix. 

When  the  number  of  points  allowed  in  a  route  exceeds  two,  some 
interesting  questions  arise  a:,  to  exactly  how  the  algorithm  should  be 
implemented.  As  an  illustra  ion,  suppose  that  in  the  e:. ample  we  allow  a 
vehicle  to  delive-  to  at  most  three  points  rather  than  two.  Now  suppose 
that  we  start  with  the  solution  J ^  in  Table  2.  We  see  that  the  two-at-a- 
t ime  potential  savings  are  all  non- pos  i  t i ve.  However,  if  we  ignore  this 
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fact  and  proceed  with  the  algorithm,  we  would  put  points  1  and  4  in  the 
same  route  since  their  potential  savings  of  zero  is  maximum.  If  we  then 
consider  iddin>  a  third  point  to  this  route,  we  find  that  adding  point  6 
has  a  potential  savings  of  2  whicli  is  maximum.  Finally,  we  combine  points 
2  and  i  i >to  a  single  route  since  this  has  a  potential  savings  oi  zero. 

The  resulting  routes  (4,15)  and  (2,3)  has  a  length  of  10.  Therefore, 
it  appears  that  we  can  get  better  information  b'-’  recalculat  ing  the  poten¬ 
tial  savings  after  each  augmentation  of  a  route.  This  recalculation  is 
not  done  in  most  implementations  of  the  C-W  algorithm. 

Clearly,  vhen  constructing  a  route  containing  more  than  two  points 
one  must  decide  where  in  the  route  to  put  each  additional  point.  The  poten¬ 
tial  savings  can  be  determined  exactly  only  by  solving  a  travelling  sales¬ 
men  problem  over  each  new  point  which  is  a  candidate  to  be  added  to  the 
route.  This  is  computationally  expensive  if  a  route  can  contain  a  large 
number  of  points.  In  most  implementations  of  the  C-W  algorithm,  new  points 
are  simpl  adde  i  on  to  the  end  of  the  route  being  constructed.  When  the 
algorithm  is  implemented  interactively  using  computer  graphics,  it  appears 
that  the  human  can  perform  an  important  role  both  in  selecting  candidate 
points  and  in  inserting  them  logicaLly  into  routes. 

Location  -  Allocation  with  Pricing 

Another  procedure  which  Krolak  and  Nelson  [5]  have  found  effective 
in  approaching  ; he  delivery  problem  utilizes  the  location  -  allocation 
model  of  Cooper  [4],  This  basic  concept  can  also  be  used  in  conjunction 
with  Theorem  1  to  generate  intuitively  appealing  columns  to  add  to  the 
set  partitioning  model. 

To  illustrate  the  procedure  consider  the  deliver}  example  illustrated 
in  Figure  2.  Again  the  circles  represent  delivery  points  and  the  square 
represents  the  depot.  The  basic  idea  is  to  use  a  suriogate  distance 
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rather  than  the  actual  distance  in  determining  the  points  which  are  assigned 
to  each  vehicle.  The  surrogate  distance  is  obtained  by  assuming  that  the 
vehicle  trivels  from  the  depot  to  a  specified  cluster  point,  represented 
in  Figure  l  by  the  dashed  circles.  It  then  makes  the  deliveries,  returning 
after  each  delivery  to  the  cluster  point.  After  all  deliveries  have  been 
made,  the  vehicle  returns  to  the  depot.  Under  this  surrogate  distance, 
the  problem  becomes  one  of  locating  the  cluster  points,  one  for  each 
vehicle,  and  then  assigning  the  delivery  points  to  each  vehicle. 


If  (a^.b^)  represents  the  coordinates  of  the  depot,  (a^,b^)  represents 
the  coordinates  of  delivery  point  i,  (x.,y.)  represents  the  coordinates 
of  cluster  point  j,  and  assuming  Euclidean  distance,  the  problem  can  be 


modeled  as  follows: 


m  n  .  .  j 

mln  l  £2((x.-a .)  +  (y.-b  )  ]^z 

x,y,z  i=0  j=l  3  J 


ij 


III 

s-t.  I  z  <  K  j  =  1,2 . m 

i=l 

n 

7  z  =  1  i=J,2,...,m 

J-l  1J 

z . .  =  0  or  1  all  i , j 

ij 

where  n  i  ■:  the  number  of  vehicles  and  K  is  the  vehicle  capacity.  When 
Zi)  "  ^ ’  lelivery  point  i  is  assigned  to  vehicle  j  and  when  z  =  0, 
delivery  point  i  is  not  assigned  to  vehicle  j. 

While  there  is  no  method  for  efficiently  solving  this  problem 
optimally,  the  following  is  an  attractive  heuristic.  Pick  a  set  of 


Figure  2.  Delivery  Example  to  Illustrate  the 
Location  -  Allocation  Model  for  Clustering. 
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locations  for  the  cluster  points.  With  these  coordinates  fixed,  solve 
the  resulting  assignment  problem.  With  these  values  of  fixed, 

solve  the  resulting  location  problem.  Continue  alternating  between 
the  assignment  and  location  problems  for  some  specified  number  of 
iteration;  or  until  there  is  no  further  improvement  in  the  objective. 
Once  a  cluster  has  been  determined,  the  vehicle  is  then  routed  among 
the  points  of  the  cluster. 

A  slight  modification  of  this  model,  together  with  Theorem  1, 

allows  us  to  generate  attractive  new  columns  for  the  set  partitioning 

problem.  Suppose  that  we  have  a  solution  to  the  set  partitioning 

problem  and  a  set  of  row  prices  P, ,p~,  .  •  .  ,P  •  Now  consider  the  model 

1  z  m 


min 

x,y,z 


m  ii  ^  o  U 

l  l  2[(x  -a  )  +  (y.-b  )  ]  z. . 

i=0  j=l  1  3  J 


“  I 


n 

l  p 


i-l  J-l 


lZij 


m 

s.t.  I  Z  <  K  j  =  1,2 . n 

i-l  3 


n 

y  z  <  1  i  =  1,2, .  .. ,m 

j-i  11 


z  =  0  or  1. 
ij 


Any  cluster  generated  by  this  model  will  have  a  positive  potential 
savings,  with  respect  to  the  surrogate  distance.  Hence,  it  corresponds 
to  an  attractive  column  to  add  to  the  set  partitioning  problem  (consider¬ 
ing  surrogate  distances). 

Since  the  second  constraint  has  been  changed  to  an  inquality,  not 
all  delivery  points  will  be  assigned  to  cluster  points.  This  simply 
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means  that  the  current  row  price  for  the  delivery  point  is  more  attrac¬ 
tive  than  the  cost  of  serving  the  delivery  point  in  alternative  clusters 
considered  by  the  model.  The  model  can  be  solved  using  the  same  heuristic 


discussed  for  the  earlier  location  -  allocation  model. 

Dlal-A-Riti e  Problem 

The  dial-a-ride  problem  is  a  specialization  of  the  general  pick-up 
and  delivery  problem.  It  is  a  much  more  complex  routing  problem  than  the 
simple  delivery  problem.  In  the  dial-a-ride  problem  we  are  given 
an  origin-destination  trip  matrix  and  an  underlying  network  on  which 
the  trips  are  tc  be  mede.  There  is  a  single  item  (people,  goods, 
etc.)  (demand  for  service)  at  each  origin  that  needs  to  be  transported 
to  its  specified  destination.  The  items  are  transported  from  origins 
to  destinations  on  vehicles  each  having  capacity  K.  We  wish  to  satisfy 
the  trip  requirements  while  travelling  the  minimum  distance. 

This  is  a  "static"  version  of  the  dlai-a-riie  problem  since  time 
is  not  considered.  There  are  a  number  of  more  complex  versions  of 
this  i  robletn,  but  this  version  is  sufficient  to  demonstrate  the  basic  ideas 
of  our  approach. 

The  set  partitioning  model  for  the  dial-a-ride  problem  is  analogus 
to  that  of  the  delivery  problem,  but  here  rows  represent  trips  rather 
than  simpl;  delivery  points.  The  vehicle  capacity  constraints  are  handled 
by  generating  only  routes  which  satisfy  them. 

Oecomposition 

In  delivery  problems,  representative  of  the  "one-i nded"  class  of 
routing  and  scheduling  problems,  we  need  only  be  concerned  with  a  single 
stop  for  a  vehicle  to  satisfy  a  particular  demand  for  service.  In  con¬ 
trast,  the  "two-ended"  class,  which  Includes  dial-a-ride  problems,  requires 
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two  stops  in  a  specific  order  (sequence)  to  satisfy  a  specific  demand  for 
service.  It  is  this  requirement  for  sequencing  of  "pickup"  and  "dropoff" 
point  pairs  which  adds  greatly  to  the  difficulty  of  handling  the  two-ended 
class  of  vehicle  routing  problems. 

When  one  considers  examples  of  two-ended  problems,  it  becomes  imme¬ 
diately  apparent  that,  the  sequencing  re<;u  irements  griatly  inhibit  the 
complex  pattern  processing  abilities  of  the  human.  Figure  1  illustrates 
an  example  of  a  25  trip  dial-a-ride  problem.  Rather  than  displaying 
order  and  structure,  the  problem  resembles  so  much  spaghetti.  It  is 
clear  that  for  such  problems  the  human  interactor  needs  more  help  in 
generating  good  columns  for  the  set  partitioning  modeL. 

Unfortunately,  ;t  is  also  more  difficult  to  apply  straight  forward 
methods  such  as  the  savings  approach  discussed  earlier  for  the  delivery 
example.  In  generating  a  dial-a-ride  route  one  must  be  concerned  with 
where  botu  the  origin  and  the  destination  occur  in  the  sequence  in  order 
to  calculate  the  potential  savings.  In  addition,  the  capacity  constraint 
may  negate  what  otherwise  appears  to  be  good  positions  for  the  origin 
and  destination  in  the  sequence. 

Because  of  this  complexity,  it  is  helpful  to  "decompose" 
the  problem  into  two  Levels  which  we  call  "clustering"  and  "chaining  . 

In  essence,  we  :onsider  a  route  to  be  made  up  of  two  components.  Clusters 
correspond  to  trips  which  can  all  be  on  a  vehicle  at  one  time,  while 
chains  correspond  to  movement  from  the  end  of  one  cluster  to  the  beginning 
of  t lie  next.  Figure  4  provides  an  example  of  five  good  clusters,  while 
Figure  5  illustrates  one  way  in  which  these  five  clusters  might  be  linked 
into  two  chains. 

The  partitioning  model  and  pricing  concepts  can  be  effectively  exploited 
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Figure  3.  An  Example  of  a  25  Trip  Dial-a-Ride  Problem 
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to  aid  in  generating  improving  clusters  and  chains  in  a  column  generation 
approach  to  solving  two-ended  vehicle  routing  problems.  One  partitioning 
model  can  be  utilized  in  generating  improving  clusters  while  another  parti¬ 
tioning  mi  del  helps  identify  better  chains. 

In  the  partitioning  model  for  clustering,  the  columns  represent  clusters 
and  the  rows  represent  the  individual  trips  (demand  for  service).  We  shall 
demonstrate,  in  the  next  section,  just  how  the  pricing  information  from  the 
partitioning  problem  can  be  used  to  generate  additional  clusters. 

In  the  partitioning  model  for  chaining,  the  columns  represent  chains 
and  the  rows  represent  the  individual  trips.  The  function  of  the  chaining 
partitioning  model  is  to  combine  the  clusters  into  good  vehicle  routes. 

The  overall  solution  procedure  consists  of  linking  the  clustering 
models  and  chaining  models  together  in  an  interactive  manner.  The 
clustering  models,  partitioning  matrix  and  pricing  information  are  used 
to  generate  good  clusters.  These  clusters  are  then  passed  to  the  chain¬ 
ing  phase  where  they  are  linked  together  via  the  chaining  models,  par¬ 
titioning  matrix  and  associated  pricing  information.  After  chaining, 
it  is  possible  to  return  to  the  clustering  phase  to  identify  additional 
clusters . 

The  next  tvi  sections  discuss  the  specifics  of  clustering  and  chaining 

Clustering 

In  the  previous  section  we  introduced  the  clustering  concept.  In 
this  section  we  shall  provide  more  details  of  the  concept  as  well  as 
the  structure  and  operationof  various  clustering  models.  Figure  6 
depicts  a  typical  cluster  (in  this  case,  three  trips). 

Clustering  makes  sense  if  the  origins  are  reasonably  close  to¬ 
gether  and  the  destinations  are  also  close  together.  One  way  to 
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develop  an  evaluation  of  such  c lrcumstanc e  is  to  locate  the  centroid 
of  the  origins,  the  centroid  of  the  destinations  and  accumulate  the 
resulting  distances  from  the  original  trips.  'a  Figure  b  we  could 
evaluate  the  distances  represented  by  (a+b+v)  +  (d+e+f).  If  this 
sum  is  small  then  it  would  make  sense  to  cluster  the  trips. 

By  utilizing  surrogate  distances  we  lose  the  actual  route  dis¬ 
tance  evaluation;  however,  we  gain  the  ability  to  evaluate  large 
numbers  of  cluster  possibilities  conveniently  and  simultaneously. 

In  Figure  6  we  might  employ  Euclidean  distances.  In  this  case  we 
could  compare  the  sum  of  the  row  prices,  p^  +  +  ^2'  8eneratec*  in 

the  covering  model  for  clustering  to  the  quantity  2(a+b+c)  +  2(d+e+f) 

+  g  to  determine  whether  clustering  is  appropriate.  We  can  think  of 
the  latter  quantity  as  a  surrogate  for  the  vehicle  routing  distance. 

We  can  develop  a  straight  forward  extension  of  the  Location  - 
Allocation  model  discussed  in  Section  7  to  identify  good  clusters 
for  the  dial-a-rlde  problem.  In  place  of  the  quantity 

IpHy./Myj-vV 

we  would  einpLoy  a  quantity 


+  l  7  2[<xj-a1)2  +  ztJ 


where  x,,  v.,  a.  and  b^  correspond  to  origins  and  x^ ,  y^ ,  and  a^ 
bj  correspond  to  destinations  of  trips. 

It  is  also  possible  to  develop  a  savings  approach  to  clustering 


in  a  similiar  manner  to  that  described  earlier  for  the  delivery  problem. 
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Chaining 

Upon  termination  of  the  clustering  process,  a  number  of  reasonably 
good  clusters  are  available.  This  set  includes  not  only  the  "best" 
cluster  sets  as  selected  by  the  covering  model  for  clustering,  but 
also  a  number  of  other  clusters  which  might  be  nearly  as  good,  but 
which  were  not  selected  in  the  optimal  solution  to  the  covering  model. 
Figure  7  illustrates  a  set  of  clusters  which  might  result  from  the 
clustering  process.  The  clusters  in  the  figure  represent  only  the 
optimal  solution  to  the  clustering  process.  In  addition,  we  might 
have  other  clusters  available,  e.g.,  a  cluster  containing  only  trips 
1  and  2. 

The  clusters  obtained  represent  good  possibilities  for  segments 
(legs)  of  a  vehicle  route.  The  next  step  in  the  process  is  to  link 
("chain")  these  clusters  into  complete  vehicles  routes.  Figure  8 
illustrates  the  chaining  concept.  Trips  1  and  2  form  one  cluster, 
while  trips  3,  4  and  5  form  another  cluster  (see  Figure  8a).  In 
Figure  8b,  trips  1  and  2  are  replaced  by  a  single  pseudo  (cluster) 
trip,  as  is  also  the  case  for  trips  3,  4  and  5.  These  cluster  trips 
are  then  chained  together. 

The  interpretation  of  chaining  is  that  a  single  vehicle  will 
service  the  first  set  of  trips  (in  Figure  8b  these  would  be  trips  1 
and  2)  and  then  proceed  to  service  the  next  set  of  trips  (i.e.,  trips 
3,  4  and  5)  in  the  chain.  Figure  9  illustrates  the  likely  vehicle 
route  to  service  the  two  clusters. 

We  shall  demonstrate  how  mathematical  models  can  be  developed  to 
aid  the  human  in  developing  good  chains  (vehicle  routes). 


Figure  7.  Example  of  Clusters  Resulting  from  the  Clustering  Process 
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a.  Two  Typical  Clusters 


b.  The  Associated  Cluster  Arcs  Chained  Together 


Figure  8.  An  Example  of  Chaining  Clusters  Together 


i 


i 

i 


!< 


52 


The  most,  fundamental  model  of  cliaining  is  based  on  network  flows. 

In  a  network  we  represent  each  cluster  trip  by  two  nodes  (an  origin 
and  1  destination  node)  and  an  arc.  We  then  add  additional  arcs  to  the 
network  which  link  tin  destination  node  of  one  cluster  arc  to  the  origin 
node  of  another  cluster  arc  on  the  basis  of  proximity  considerations. 

In  Figure  10,  we  indicate  one  network  flow  model  for  the  clusters  of 
Figure  7. 

To  establish  the  flow  variables  for  the  network  flow  problem  consider 
the  case  where  each  trip  appears  in  exactly  one  cluster.  We  assign  every 
arc  an  upper  capacity  of  1,  the  dashed  arcs  a  lower  capacity  of  0  (zero) 
and  the  solid  cluster  arcs  a  lower  capacity  of  1  indicating  that  these 
trips  must  be  service!.  Arc  flow  costs  for  the  network  flow  problem 
•ire  the  associated  vehicle  travel  distances.  If  the  "starting  node' 
and  the  "ending  node"  represents  the  vehicle  storage  depot,  then  we  may 
ais(  assi>  n  additional  costs  to  the  arcs  originating  at  the  starting 
node  to  reflect  the  fixed  cost  of  using  each  vehicle  (provided  they 
are  all  the  same). 

If  the  network  flow  problem  for  chaining  contains  no  circuits  and 
if  the  sane  trip  does  not  appear  in  two  different  clusters  in  the  network, 
then  the  resulting  chains  (vehicles  routes)  will  be  valid  ones.  If  either 
of  these  two  conditions  are  explicitly  included  in  the  model,  then  the 
underlying  network  flow  problem  becomes  considerably  more  complicated. 

We  can  circumvent  some  of  the  problems  associated  with  the  difficult 
problems  of  chaining  by  utilizing  a  partitioning  model  of  chaining.  (Actually 
we  substitute  one  difficulty  for  another.)  In  a  partitioning  model  of  chaining. 
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wc  issori  ltf  a  column  of  tin-  partitioning  matrix  lot  each  It-asible  chain 
(vt'  ilde  route)  and  a  row  for  each  trip.  Column  j  contains  1  in  row 
i  i:  trip  i  is  serviced  hy  chain  (vehicle  route)  j.  Otherwise,  it 
contains  a  0  (zero) .  The  zero-one  variable  associated  with  the 
columns  provide  indications  of  which  chains  were  selected  in  the  optimal 
partitioning  solution.  Table  1  illustrates  a  partitioning  matrix 
tor  several  chains  in  Figure  LO. 

besides  solving  the  partitioning  model  optimally  (which  we  would  not 
likely  do),  the  obvious  difficulty  with  such  a  model  is  in  generating 
the  candidate  chains  (columns).  To  aide  in  this  column  generation 
process,  we  could  employ  (1)  the  human  (which  would  certainly  be  part 
of  any  interactive  process),  (2)  a  savings  approach  to  combine  clusters 
into  chains,  (3)  a  network  flow  model  or  (4)  some  other  method.  The 
first  two  if  these  column  generation  approaches  should  be  reasonably 
intuiiive  to  develop.  We  shall  briefly  describe  the  third  one. 

Suppi'se  we  have  identified  some  candidate  columns  (chains)  for  the 
partitioning  problem  (e.g.  the  individual  clusters  f roia  the  clustering 
process).  Then,  the  partitioning  model  will  yield  (1)  an  optimal 
solution  and  (2)  a  set  of  row  prices.  As  before,  a  column  j  not  in  the 
solution  appears  favorable  it 


P  .  a  .  . 

ill 


c.  0. 
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N'i'w,  r  is  the  i-nsl  (il  i  s  t  a  nr  r )  ol  servicing  the  chain.  The  )  p.a.. 

i  .  J  t  I 

is  I  lie  ( -list  ol  servicing  the  given  trips  in  the  current  partitioning 
solu!  ion.  Consider  tiow  these  two  terms  might  he  represented  in  the 


k  of 

!■'  i  gu  i 

,  Hi. 

upposc  wo 

associate  with  each  dashed 

arc 

the 

neg 

1 1  i  vc 

ol  what- 

lists 

that 

ire  incurred  in  traversing 

the 

arc 

(th 

is  is 

t he  same  as 

before  with  a  sign  change).  Associated  with  each  clustering  arc 
the  quant i tv  1  p  -  c ,  where  the  term  J  p  represents  the  sum  of  row 
prices  for  trips  in  the  cluster  and  c  represents  the  cost  of  the  cluster. 
■  ith  these  costs  defined  on  the  network,  we  seek  a  pain  (or  paths)  in 
the  network,  from  the  starting  to  the  ending  node,  which  has  positive 
total  cost.  It  should  he  clear  that  such  a  path  satisfies  the  condition 
lor  a  potentially  improving  chain,  and  can  be  added  to  the  partitioning 
prob l em . 

Again,  with  this  network  flow  model  (a  shortest  path  model  for  a 
sinrlc  chair)  we  have  tie  difficulties  associated  with  circuits  and 
with  t  lie  same  tri,  serviced  several  times  by  a  single  chain.  The 
phenomenon  of  a  trip  being  serviced  several  times  seems  to  occur  in¬ 
frequently  since  in  seeling  a  least  cost  solution,  the  model  tends  to 
avoid  such  an  instance.  When  it  does  occur,  an  attractive  heuristic 
is  to  simply  delete  one  of  the  conilicting  clusters  from  the  column 
being  generated.  We  are  currently  utilizing  the  human  interactor  to 
handle  the  case  where  the  flow  problem  contains  circuits.  The  human 
breaks  the  circuit;  and  patches  the  paths  back  together.  If  we  re- 


57 


■a  riot  In  pi, 'cedine  to  locating  a  si  little  improvi  ng  path  each  time, 

(  •  ■  if,  1 1 1  i  tip  shortest  path  model  lends  itself  more  intuit  ively  to 
[  ,e  Jev  I  opiiu  it  of  good  heuristic  procedures  for  solving,  hoth  ot  the 
■•i,  or  i  M  in  Hies  associated  with  tlu  flow  mode.  in  particular  it 
i>  mile n  easier  tor  the  hu nan  to  break  the  circuits  in  a  single  path 
than  in  ;nu  1  t  i  ;>a  t  hs . 

Conclusions 

An  interactive  el  ia  1-  l  -r  id.;  system  based  on  the  concepts  presented 
here'  has  been  implemented  on  a  Chromatics  color  graphics  terminal  inter¬ 
laced  w'lh  a  CYBKR  74  mainframe  computer.  Although  the  system  is  in 
many  wavs  ver /  rudimentary,  it  dramatically  indicates  the  potential 
lor  this  kind  of  interactive  optimization.  We  are  in  the  process  of 
making  extensive  modifications  in  the  software  and  are  testing  a 
variety  of  now  models  and  heuristics  to  aid  in  route  generation. 
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IV.  ORUI-R  PICKING  I'KOBLKM 

One  >>f  tin'  nost  l  und.imfiita  I  ,  among  r  li<*  man;.  impu  r  Gin  t  i  Lem  tclricv- 
a  I  problem?,  associated  with  warehousing  and  mater  ial  •.  dandling  j  >.  what 
we  will  call  the  simple  'order-picking  problem."  An  order  consists  of  a 
subset  of  the  items  si  >red  in  a  warehouse.  When  an  order  is  requested, 
a  vehicle  is  dispatched  from  the  shipping  area  to  collect  or  "pick"  the 
items  in  the  order  and  transport  them  back  to  the  shipping  area.  The 
objective  is  to  minimize  distance  traveled  by  the  vehicle.  We  will 
assume  that  a  vehicle  picks  only  one  order  at  a  time  and  that  an  order 
does  not  exceed  the  vehicle's  capacity. 

For  an  arbitrary  aisle  configuration  within  a  warehouse,  the  order¬ 
picking  problem  is  easily  recognized  as  a  variant  of  the  well  known  and 
difficult  t  :>  sol' e  traveling  salesman  problem.  Fortunately,  the  most 
common  warenou.se  aisle  configuration  is  that  given  in  Figure  1.  For 
this  configuration  we  will  present  an  algorithm  for  the  order-picking 
problem  which  is  linear  in  the  number  of  aisles. 

Graph  Kepjm  sen  Ur  on 

Consider  an  order  containing  in  items  which  is  to  be  picked  in  a 
warehouse  with  the  aisle  configuration  illustrated  in  Figure  1.  De¬ 
fine  a  graph  C  by  associating  a  vertex  v^  with  the  shipping  area  location, 
a  vertex  with  the  location  of  each  item  i  =  1,  2,,,.,tn  in  the  order, 
and  vertices  a.  and  b,  with  the  ends  of  each  aisle  j  =  1,  2,...,n. 

Connect  any  two  vertices  in  G  which  correspond  to  adjacent  locations  in 
the  warehouse  by  an  unlimited  number  of  parallel  arcs,  each  with  length 
equal  to  the  direct  distance  between  the  two  locations.  An  example 
graph  is  given  in  Figure  2  corresponding  to  the  warehouse  and  order 


Figure  2:  Order-Picking  Graph  Where  Each  Arc  Represents 
An  Unlimited  Number  of  Parallel  Arcs 
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depicted  in  Figure  1.  For  simplicity  the  parallel  arcs  connecting 
adjacent  vertices  have  been  represented  by  a  single  arc.  Also  for 
simplicity  v^  is  depicted  as  being  between  vertices  and  in¬ 
stead  of  coinciding  with  B. . 

4 

An  order  picking  tour  is  a  cycle  in  G  which  includes  each  of  the 

vertices  v.  for  i  =  0.  at  least  once.  Note  that  since  each 

i 

adjacent  pair  of  vertices  is  connected  by  an  unlimited  number  of  paral¬ 
lel  arcs,  we  can  assume  without  loss  of  generality  that  an  order-picking 
tour  contains  any  arc  at  most  once.  The  order-picking  problem  is  then 
to  find  a  minimum  length  order-picking  tour  in  G. 

A  subgraph  Tc  G wi 11  be  called  a  tour  subgraph  if  there  is  an  order¬ 
picking  tour  which  uses  each  arc  in  T  exactly  once.  The  following  cha¬ 
racterization  of  a  tour  subgraph  is  a  specialization  of  a  well  known 
theorem  on  Euler  graphs. 

Theorem  1.  A  subgraph  TcG  is  a  tour  subgraph  if  and  only  if  (a)  T  con¬ 
tains  all  vertices  v^  for  1  =  0,  1,  2,.,.,m;(b)  T  is  connected;  and 
(c)  every  vertex  in  T  has  even  degree. 

Given  a  tour  subgraph,  we  will  show  in  a  later  section  that  an 
order-picking  tour  can  be  very  efficiently  determined.  Hence,  by 
developing  an  efficient  procedure  for  finding  a  minimum  length  tour 
subgraph  we  can  efficiently  solve  the  order-picking  problem.  The 
following  corro . laries  of  Theorem  1  are  useful  in  characterizing  a 
tour  subgraph. 

Corrollary  1,  A  minimum  length  tour  subgraph  contains  no  more  than 
two  arcs  between  any  pair  of  vertices. 

Corrollary  2.  If  (P,P)  is  any  node  partition  of  a  tour  subgraph,  there 
is  an  even  number  of  arcs  with  one  end  in  P  and  the  other  in  P. 


let  Lj  be  a  subgraph  of  G  consisting  of  vertices  and  and 

everything  in  C  to  the  left  of  a.  and  b  but  not  including  arcs  and 

vertices  ietweeu  or  to  the  right  ot  a.  and  b,.  An  illustration  is 

J  J 

given  in  l  igure  S.  A  subgraph  T.iJLj  in  an  L^-partiaL  tour  subgraph 

if  it  contains  every  v .  e  Lj  and  at  least  one  of  {a^,b^},  every  vertex 

in  T.  except  possible  for  a.  and  b.  have  even  degree,  and  T.  has  either 

J  1  J  J 

a  single  connected  component  or  two  connected  components  with  in 
one  component  and  in  the  other.  Some  examples  illustrating  L^-partial 
tour  subgraphs  are  given  in  Figure  4. 

1  2 

Two  Lj-partial  tour  subgraphs  and  are  equivalent  if  has 

the  same  degree  in  both,  bj  has  the  same  degree  in  both,  and  either  both 
have  one  connected  component  or  both  have  two  connected  components. 

1  2 

Theorem  2.  If  T.  and  T.  are  two  equivalent  L.  partial  tour  subgraphs, 

- - -  J  J  J 

and  R  is  a  subgraph  of  G-L .  such  that  T^  U  R  is  a  tour  subgraph,  then 

2 

Tj  U  R  is  also  a  tour  subgraph. 

Proof.  Since  T . *  II R  is  a  tour  subgraph,  R  must  contain  all  v.  c  G-L. 

-  1  or.  i  j 

and  except  possibly  for  and  b  ,  each  vertex  in  R  must  have  even  degree. 
2 

Hence,  f.  y  K contains  all  v.  t  G.  Since  a,  and  b.  have  the  same  de- 
J  J  J  J 

2  2 

gree  in  both  T.  and  T.  and  all  other  vertices  in  T.  are  even,  all 
l  J  J 

> 

vertices  in  T ^ “  'U  R  must  hive  even  degree. 

Now  note  that  shrinking  a  connected  component  of  a  graph  to  a  single 


vertex  does  not  affect  the  connectivity  of  the  graph.  If  in  T.  (J  F.  and 

2  1  2  J 
Tj  y  R  we  shrink  the  connected  components  of  T^  and  T^  to  single  vertices 

1  2 

the  resulting  graphs  are  identical  (i.e.,  either  Tj  and  T^  each  shrink 

to  a  single  vertex  connected  to  a.  and  b,  in  the  same  fashion  or  they  each 

J  j 

shrink  to  two  vertices  connected  to  a.  and  b,  in  the  same  fashion).  There- 

J  j 

1  2 

fore  since  T.  II  R  Ls  connected  T.  II  R  is  also  connected.  It  then  follows 
J  j  u 

2 

from  Theorem  1  that  Tj  is  a  tour  subgraph. 


Cons  t  rue  t  i  ng  Minimum  Length  lour  Subgraphs 

Before  constructing  a  minimum  length  tour  subgraph,  we  can 
simplify  the  order-picking  graph  somewhat  by  use  of  the  following 
resu i t . 

Lemma  1:  If  an  aisle  does  not  contain  at  least  one  item  in  a  given 
order  then  there  is  an  optimum  order-picking  tour  which  does  not 
traverse  the  aisle. 

Proof:  Consider  anv  two  vertices  v.  and  v.  which  are  adiacent  in  an 

-  !  j 

optimum  tour  (see  Figure  5).  The  minimum  distance  traveled  between 

v.  and  v.  is  D  -  min  {a+c,b+d}  +  h.  Clearly,  this  can  be  attained 
l  J 

without  traversing  any  aisle  other  than  those  containing  v^  and  v^ . 

Therefore  we  can  delete  from  consideration  those  arcs  corresponding 
to  aisles  which  contain  no  items  in  the  order.  For  simplicity,  in 
the  following  discussion  we  will  assume  that  this  has  been  done  and 
the  aisles  and  corresponding  vertices  renumbered. 

Consider  any  minimum  length  tour  subgraph  T  and  the  vertex  pair 

a.,b  corresponding  to  the  ends  of  some  aisle  j.  From  Clorrollarv 

r  ! 

]  we  know  that  the  degree  of  everv  vertex  in  T  is  even  and  from  Cor— 


rollarv  2  we  know  that  for  anv  j  >  1  the  number  of  arcs  connecting 

a.  .  and  a  plus  the  number  of  arcs  connecting  b,  ,  and  b.  is  also 
J-1  j  j-1  J 

even.  Eliminating  all  possibilities  which  do  not  satisfy  these  condi¬ 
tions,  we  are  left  with  the  80  cases  in  Figure  6  for  the  connectivity 

of  a  and  in  T.  Vertices  between  a.  and  b,  have  been  ignored  except 
J  1  J  J 

where  they  are  necessary  to  indicate  connectivity. 

From  Theorem  1  the  degree  of  each  vertex  in  T  must  be  even.  This 


) 


Figure  6:  Continued 
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order.  Hi-na;,  we  can  eliminate  eases  D.  through  I)  since  they  would 

I  6 

nev-r  appear  in  an  optimum  tour.  For  example,  case  A(^  dominates  cases 

0  .  D  ,  and  Dc  since  for  those  cases  we  do  not  need  Lo  visit  a.. 

J  1  >  J 

Similarly  case  A.,  dominates  cases  D„  ,  D.  ,  and  D,. 

1  /  J  4  6 

Note  that  between  vertices  a.  and  b.  the  tour  subgraph  takes  one 

J  J 

of  the  forms  in  Figure  7.  For  (a),  (b),  (c) ,  and  (e)  the  distance 
traveled  within  the  aisle  is  unique.  For  an  aisle  containing  p  items 
to  be  picked,  there  are  p-1  possibilities  corresponding  to  case  (d) . 
However,  the  only  possibility  that  we  need  to  consider  is  determined  by 
the  two  adjacent  items  within  the  aisle  which  are  fartherest  apart.  As¬ 
suming  p  items  to  be  picked  in  the  aisle,  this  requires  p-1  comparisons. 
We  will  assume  that  this  has  been  done  for  each  aisle.  Hence,  the  length 
of  that  part  of  the  tour  between  and  Ik  is  known  for  each  of  the  arcs 
A1  through  A^Q. 

From  Figure  6  we  see  that  for  i  =  2,3,...,n-l  there  are  only 
four  degree  combinations  for  and  b^  in  a  L^-partial  tour  subgraph 
which  can  lead  to  a  tour  subgraph.  Denoting  these  combinations  by 
(degree  a.,  degreea  )  they  are  (1,1),  (0,2),  (2,0),  and  (2,2).  If  we 
examine  the  connectivity  of  each  of  these  combinations,  it  follows  from 
the  definition  that  the  (0,2)  and  (2,0)  L^-partial  tour  subgraphs  must 
have  only  one  connected  component.  The  (1,1)  L^-partial  tour  subgraphs 
must  also  have  only  one  connected  component  since  to  have  two  components 
with  one  connected  to  a^  and  the  other  to  b^  implies  that  some  vertex 
other  than  a.  and  b^  has  odd  degree.  Finally,  the  (2,2)  ^-partial  tour 
subgraphs  can  have  either  one  or  two  components.  Hence,  there  are  five 
equivalence  classes  which  we  will  denote  (1,1),  (0,2)  (2,0),  (2, 2, 1C), 
and  (2,2, 2C).  Under  our  definition  of  'equivalent/'  all  possible  Lj- 


) 
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partial  tour  subgraphs  fall  into  one  of  these  classes. 

It  follows  from  Theorem  2  that  we  only  need  to  concern  ourselves  with 
the  minimum  length  L^-partial  tour  subgraph  in  each  equivalence  class  in 
constructing  a  minimum  length  tour  subgraph.  To  see  how  a  minimum  length 
Lj-partial  tour  subgraph  can  be  constructed  for  each  equivalence  class, 
assume  that  for  some  j  =  3,4,...,n-l  we  have  minimum  length  L .  ^-partial 
tour  subgraphs  for  each  of  the  five  equivalence  classes. 

First  consider  the  minimum  length  L.  ^-partial  tour  subgraph  corres¬ 
ponding  to  the  (1,1)  class.  If  we  add  the  arcs  indicated  in  Figure  6  case  A^ 

between  a.  ,  and  b.  ,,  between  a.  ,  and  a.  and  between  a.  .  and  b.,  we 
J-l  J-1  I"1  2  '  J"1  1 

will  have  constructed  an  L^-partial  tour  subgraph  in  the  (2, 2, 1C)  equivalence 
class.  Figure  7  indicates,  in  terms  of  the  cases  enumented  in  Figure  6 
all  of  the  possible  transitions  from  a  L..  ^-partial  tour  subgraph  to  a 
Lj-partial  tour  subgraph.  Note  that  the  cost  of  the  transition  is  simply 
the  length  of  the  arcs  added  in  going  from  an  j-1  equivalence  class  to  a 
j  equivalence  class.  For  situations  such  as  the  transition  from  (1,1)  to 
(1,1)  in  Figure  7  the  cost  of  the  transition  is  the  minimum  of  the  cases 
A^,  A^,  A^.  Actually  the  A^  case  is  dominated  here  as  are  A  , 

and  A^q  for  the  (2, 2, 1C)  to  (0,2),  (2,2,1C)  to  (2,0),  and  (2, 2, 1C)  to 
(2, 2, 1C)  transitions  respectively.  Hence,  these  cases  can  be  deleted  from 


consideration. 

Therefore,  to  find  the  minimum  length  L^-partial  tour  subgraph  for 
each  equivalence  class,  given  the  minimum  length  I,  ^  partial  tour  subgraph, 
we  first  add  the  Length  of  the  partial  tour  subgraph  for  each  class 

to  the  transition  cost.  Then  for  each  j  equivalence  class,  we  pick  the 
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minimum.  the  sums  which  correspond  to  transitions  to  that  class. 

For  ]  -  1  and  j  =  n  the  possible  cases  are  shown  in  Figure  8. 

Aisle  1  has  only  a  single  equivalence  class.  The  transition  is  shown 
in  Figure  9.  The  cost  of  a  minimum  length  I^-partial  tour  subgraph 
for  each  equivalence  class,  is  simply  the  transition  cost.  lienee  we 
can  easily  start  the  process.  Once  a  minimum  length  L^-partial  tour 
subgraph  has  been  determined  for  each  equivalence  class,  the  transitions 
indicated  in  Figure  10  take  us  to  a  tour  subgraph. 

We  can  couch  this  procedure  in  dynamic  programming  terms  or 
equivalently  consider  it  a  shortest  path  problem  in  an  acylic  graph 
(i.e.  dirt  t  all  of  the  arcs  from  lower  to  higher  numbered  aisles) 
where  the  vertices  in  Figures  7,  9,  and  10  correspond  to  states. 

Tour  Construction 

We  have  developed  a  procedure  for  finding  a  minimum  length  tour 
subgraph  T  from  the  order-picking  graph.  There  remains  the  question 
of  how  to  construct  a  tour  from  T.  The  following  is  a  very  simple 
and  efficient  procedure  for  constructing  an  optimum  order-picking 
tour  f  rom  T . 

Step  1. . 

Step  2. 

Step  3. 

Step  4. 


Begin  the  tour  by  letting  vQ  be  the  first  vertex  visited. 

Let  v  be  the  vertex  currently  being  visited. 

If  there  is  a  pair  of  unused  parallel  arcs  in  T  incident  to 
v*,  use  one  of  them  to  get  to  the  next  vertex.  Go  to  Step  2, 
If  there  are  any  unused  single  arcs  in  T,  (i.e.,  not  one  of  a 
pair  of  parallel  arcs)  use  one  of  them  to  get  to  the  next 
vertex.  Go  to  Step  2. 


AISLE  n 
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Figure  10:  Transitions  From  A  Ln-Partinl  Tour 
.  Subgraph  To  A  Tour  Subgraph 


Step  5.  If  there  is  a  pair  of  parallel  ires  in  T  with  one  still  unused, 
use  it  to  get  to  the  next  vertex.  Go  to  Step  2. 

Stop  6.  Stop.  The  order-picking  tour  is  complete. 

We  first  note  that  this  procedure  constructs  a  cvcle  in  the 
tour  subgraph  T  which  starts  at  v^.  Hence,  it  must  also  stop  at 
since  each  vertex  has  an  even  number  of  edges  incident  to  it.  Every 
edge  incident  to  v^  must  be  included  in  the  cycle  or  else  the  procedure 
would  not  have  stopped.  The  only  question  left  to  resolve  is  whether 
or  not  the  cycle  includes  every  edge  in  T.  To  show  this  we  will  use  the 
following  results. 

Theorem  3;  If  ar.d  a2  are  a  pair  of  parallel  arcs  In  a  minimum  length 
tour  subgraph  T,  there  does  not  exist  two  arc  disjoint  cycles  in  T  with 
in  one  and  in  the  other. 

Proot :  Suppose  that  there  are  two  such  cycles.  If  we  delete  a^‘  and  a2  from 
T,  the  de  ;ree  of  each  vertex  in  T  is  still  even  and  T  is  still  connected. 
This  contradicts  the  assumption  that  T  is  a  minimum  length  tour  subgraph. 
Corollary  3.1:  The  tour  construction  procedure  cannot  include  in  the 
cycle  cnlv  one  of  the  pair  of  parallel  arcs  and  a2> 

Proof :  Assume  that  only  one  of  the  pair,  say  was  included  in  the.  cycle. 

We  can  then  construct  a  second  cycle,  which  includes  a2  by  starting  with 
the  vertex  at  one  end  of  a2  and  adding  incident  edges  and  vertices  not 
in  the  first  cycle  until  we  reach  the  vertex  on  the  other  end  of  From 

Theorem  3,  this,  contradicts  the  assumption  that  T  is  minimum  length  tour 
subgraph. 

Theorem  4:  Thu  tour  construction  procedure  includes  every  arc  of  T  in 


the  cycle. 


Proof  Assume  chat  some  arc  of  T  Is  not  included  in  the  cycle.  Since  I 
is  connected,  there  must  be  an  arc,  say  a,  not  included  in  the  cvelo 

which  has  at  least  one  of  its  vertices,  say  v  in  the  cycle.  The  degree 

» 

of  v  in  T  is  even  and  the  cycle  uses  an  even  number  of  arcs  incident  to 

»  » 
v  .  Hence,  there  are  an  even  number  of  arcs  incident  to  v  which  are  not 

in  the  cycle.  From  Corroilary  3.1  any  pair  of  parallel  arcs  incident  to 

v  are  either  both  in  the  cycle  or  neither  in  the  cycle.  Therefore,  if 
» 

a  is  one  of  a  parallel  pair  of  arcs,  neither  are  in  the  cycle. 

t 

Now  consider  the  last  time  the  cycle  visited  v  .  It  could  not  have 
! 

left  v  via  one  of  a  pair  of  parallel  arcs  since  this  would  violate 

f 

Corroilary  3.1  if  neither  were  in  the  cycle  already,  and  a  would  have 

taken  precedence  if  one  of  the  pair  was  already  in  the  cycle.  Hence,  it 

r 

must  have  left  via  a  single  arc.  This  implies  that  a  is  also  a  .single 

arc  since  it  woul  !  have  taken  precedence  if  it  were  one  of  an  unused  pair 

We  st e  from  Figure  6  that  a  vertex  in  T  has  at  most  two  single  arc 
! 

i n^  ide.it  to  it.  Since  v  has  incident  to  it  one  single  arc  in  the  cycle 
and  on  •  single  arc  not  in  the  cycle  and  since  the  number  of  unused  arcs 

»  I 

idcat  tc  v  i:,  even,  there  must  be  a  pair  of  arcs  incident  to  v  with 
in  !  !.<•  i  w  I  e  and  the  other  not  in  the  cycles.  From  Corroilary  3,1 
t  ::  i  .  i  ,  a  cent  rad  i.  t  ion  . 

t  <  ns  ion--,  jiinl  t  mu  lu  ions 

I  lie  procedure  presented  here  provides  a  very  efficient  procedure 
! .  ■  >!  •  orde;  piikiig  i  r  hie-.  :  i  the  a  i  .!•■  cont  igur.it  ion  given 

its  I'i  ,ure  1.  A  f  ■  v;  variations  in  :  hi  .  con;  igur.it  ion  eon  he  handled 
without  much  increase  in  computational  effort. 

The  loidirig  dock  can  he  along  .lev  of  the  aisles  without  .,  P  -i  i  i 


the  model.  It  can  also  be  at  the  end  of  the  aisles  but  not  at  an 


intersection.  This  is  most  easily  handled  by  inserting  an  'artificial 

aisle''  with  length  sufficiently  long  to  keep  it  from  being  traversed. 

The  distance  between  ^  and  does  not  have  to  be  the  same  as 

the  distance  between  b.  .  and  b,.  For  example,  the  configuration  can 

1“1  1 

be  contacted  into  a  circle  with  radial  aisles  and  still  be  amenable 
to  the  procedure. 

If  there  are  items  along  the  ends  of  aisles  as  well  as  within 

aisles,  the  basic  procedure  can  still  be  applied  but  it  is  slightly 

more  complex.  For  the  Figure  2  case,  the  number  of  arcs  in  a  tour 

subgraph  incident  to  aj_^  anc*  on  the  ri8hc  1 s  t{*e  same  as  the 

number  of  arcs  incident  to  and  b^  respectively  on  the  left.  Hence 

in  setting  up  the  states  in  Figure  7  we  only  have  states  corresponding 

to  degrees  to  the  left  of  aj_^»  hj  ay  and  there  are  items 

to  be  picked  along  the  ends  of  aisles,  the  number  of  arcs  in  a  tour 

subgraph  incident  to  a..  ^  and  ^j-l  °n  the  is  not  necessarily  the 

same  as  the  number  of  arcs  incident  to  a^.  and  b^  on  the  right.  We  can 

handle  this  by  inserting  a  set  of  states  in  Figure  7  which  correspond 

to  the  number  of  arcs  incident  to  a.,  and  b._^  on  the  right.  The 

1  l  J 

analysis  i;  then  analogous  to  that  done  for  the  previous  case.  This 
approximately  doubles  the  computational  effort. 

The  basic  ideas  can  also  be  extended  to  the  case  where  cross-overs 
are  allowed  within  the  aisles  as  well  as  at  the  ends  of  aisles.  How¬ 
ever  this  dramatically  increases  the  number  of  states  and  cases  which 
must  he  considered.  In  the  case  where  cross-overs  within  aisles  were  not 
allowed,  we  only  had  to  be  concerned  with  combinations  of  degrees  for  the 
two  vertices  cor responding  to  the  ends  of  the  aisle. 


/* 


If  we  allow  p  cross- 


H2 


overs  within  an  aisle,  then  we  have  to  consider  degree  combinations 
for  p  +  2  vertices.  In  addition,  we  only  had  to  consider  one  connectivity 
class  for  four  of  the  degree  combinations  and  two  connectivity  classes 
for  the  fifth  in  modelling  the  prohLem  with  no  cross-overs  within 
aisles.  If  we  allow  cross-overs  within  aisles  the  number  of  connectivity 
classes  also  increases.  For  example,  if  one  cross-over  is  aLlowed  with¬ 
in  each  aisle,  these  are  the  13  degree  combinations  (1,1,2),  (1,2,1), 
(2,1,1),  (2,2,0),  (2,0,2),  (0,2,2),  (2,2,2),  (0,0,2),  (0,2,0),  (2,0,0), 
(1,1,0),  (1,0,1),  (0,1,1).  The  first  six  can  have  either  one  or  two 
connected  components.  The  last  six  have  one  connected  component.  The 
combination  (2,2,2)  can  have  either  one,  two,  or  three  connected  components. 
In  addition,  to  specify  the  equivalence  class  if  there  are  two  connected 
components,  we  need  to  know  which  pair  of  the  three  intersertion  nodes 
are  in  the  same  component.  Hence  the  (2,2,2)  case  requires  five  equivalence 
classes.  This  is  a  total  of  21  equivalence  classes-  This  can  still  be 
handled  with  a  reasonable  amount  of  computational  effort.  However, 
the  procedure  does  not  seem  practical  for  more  than  two  or  three 
cros  -over  :  with  eacii  aisle. 


V.  SUMMARY 


In  our  three  areas  of  primary  focus  we  feel  that  we  have  made  a 
significant  fundamental  contribution.  We  have  developed  and  partially 
tested  interactive  methodology  for  addressing  the  fleet  scheduling  problem, 
in  addition  to  greatly  improving  t he  fleet  scheduling  process.  Insight 
gained  from  this  effort  has  provided  an  exciting  new  direction  for  research 
into  other  complex  scheduling  problems.  One  of  these  is  the  scheduling 
of  dry  bulk  carrying  ships  for  which  we  have  a  research  effort  underway. 
Another  is  the  scheduling  of  ship  construction  into  shipyards.  This  is 
an  area  that  we  are  exploring  for  research. 

Interactive  methodology  for  the  delivery  problem  has  been  developed 
and  tested.  This  methodology  seems  clearly  superior  to  previously  available 
methodology  for  addressing  a  broad  class  of  routing  and  delivery  problems. 

While  we  started  out  trying  to  develop  interactive  methodology  for 
order-picking  problems,  we  found  that  we  could  develop  methodology  to 
solve  certain  classes  of  these  problems  optimally  without  human  interaction. 
Under  the  current  ONR  contract,  we  are  exploring  other  realistic  order-picking 
problems  to  determine  if  efficient,  methods  can  also  be  developed  for  them. 

In  addition  a  number  of  new  research  areas  were  indentified  and  are 
currently  being  examined  within  the  Production  and  Distribution  Research 
Center  at  Georgia  Tech.  The  Center  was  established  under  funding  by  the 
Office  of  Naval  Research.  These  areas  include  interactive  facility  design, 
storage  and  retrieval  system  design,  and  packing  and  packaging.  The  basic 
methodology  being  developed  for  these  problems  can  be  considered  as  ex¬ 
tensions  of  interactive  concepts  and  models  generated  under  this  contract. 

In  addition  to  an  ongoing  interaction  with  the  SURFPAC  scheduler  at 
Norfolk,  we  have  been  involved  in  a  number  of  other  briefings  and  discussions 
with  Navy  personnel  related  to  problems  of  interest  to  the  Navy,  particularly 


84 


in  the  area  of  material  handling. 

On  .July  14,  1979  a  briefing  on  material  handling  systems  design  was 
provided  at  NSC-Oakland  to  approximately  36  U.  S.  Navy  personnel  from 
NSC-Oakland,  NSC-Puget  Sound,  NSC-San  Diego,  NAS-Alameda,  NARIr-Alameda , 

Mare  Island  Shipyard,  and  the  Navy  School  of  Transportation.  The  majority 
of  those  attending  were  senior  officers,  additionally,  several  supply 
officers  were  present.  The  briefing  was  provided  with  the  objective  of 
preparing  the  executive  managers  for  the  design,  development,  and  manage¬ 
ment  of  advanced  material  handling  systems  required  to  meet  the  changing 
business  requirements  of  the  U.  S.  Navy. 

On  July  15,  1979  a  meeting  was  held  with  NSC-Oakland  personnel  to 
review  ongoing  projects  involving  the  analysis  and  design  of  automated 
material  handling  systems.  The  projects  were  critiqued  and  recommendations 
were  made  for  future  analyses  and  issues  to  be  addressed . 

On  September  24-25,  1979  briefings  were  provided  in  Washington,  D.  C. 
to  NAVSUPSYSCOMHG  personnel.  The  first  briefing  was  for  the  Cliief  of 
the  Supply  Corp  and  Commander  of  the  Supply  Systems  Command  and  his 
Executive  Board,  consisting  of  Admirals  and  Captains.  Subsequently,  an 
expanded  briefing  was  given  to  NAVSUPSYSCOM  Functional  managers  and  others 
involved  in  policy-making.  Senior  level  officers  (Captains  and  Commanders) 
and  professional  level  civilians  attended  the  briefing.  Empahsis  was  placed 
on  strategic  planning,  analysis,  design,  and  evaluation. 

Researchers  who  have  worked  on  various  projects  related  to  research 
under  this  contract  include  Dr.  H.  Donald  Ratliff,  Dr.  John  A.  White,  Dr. 
John  J.  Jarvis,  Dr.  Arnie  Rosenthal  and  graduate  students  Bruce  Brownlee, 
Brian  Thorn,  Mark  Goetschalchx,  Frank  Cullen,  and  Richard  Sharp. 
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